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Abstract 



We consider the effective topological field theory on Euclidean D-strings wrapping on a 2- 
cycle in the internal space. We evaluate the vev of a suitable operator corresponding to the 
chemical potential of vortices bounded to the D-strings, and find that it reduces to the partition 
function of generalized two-dimensional Yang-Mills theory as a result of localization. We argue 
that the partition function gives a grand canonical ensemble of multi-instanton corrections for 
four-dimensional 7V=2 gauge theory in a suitable large N limit. We find two-dimensional gauge 
theories that provide the instanton partition function for four-dimensional 7V=2 theories with 
the hypermultiplets in the adjoint and fundamental representations. We also propose a partition 
function that gives the instanton contributions to four-dimensional J\f=2 quiver gauge theory. 
We discuss the relation between Nekrasov's instanton partition function and the Dijkgraaf-Vafa 
theory in terms of large N phase transitions of the generalized two-dimensional Yang-Mills 
theory. 
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1 Introduction and Summary 



One of the most significant developments in supersymmetric (SUSY) gauge theories in 
recent years is that non-perturbative dynamics including instanton corrections becomes to 
be captured by interesting physical/mathematical systems. In particular, the prepotential 
of four-dimensional J\f=2 SUSY gauge theory and the superpotential of Af=l SUSY gauge 
theory, which are both holomorphic, are obtained from the free energy of a random 
partition [1J and a random matrix model [21 El Ej , respectively. It is also pointed out that 
the random partition describing the prepotential of the M=2 theory is related to a large 
N limit of two-dimensional Yang-Mills theory [5|. Similar developments have also come 
to topological string theories. For example, amplitudes on a singular Calabi-Yau manifold 
are obtained from some statistical model like the melting crystal [HIE!- Furthermore, the 
closed topological string amplitude also has a dual open string picture which is described 
by the large N limit of three-dimensional Chern-Simons gauge theory [HJ [U] . 

An interesting fact is that these systems are "bosonic" even though the original gauge 
theories and the topological string theories are preserving SUSY. In prior to the discov- 
ery of the relationship between the holomorphic quantities of the SUSY systems and the 
bosonic partition functions, relations between SUSY gauge theories and "bosonic" inte- 
grable systems are already pointed out (see for review [TO] ) - In fact, the Seiberg-Witten 
curve of the four- dimensional J\f=2 SUSY gauge theory is especially the spectrum curve 
for the integrable system and vacua in the four- dimensional J\f=l theory are equilibrium 
points of the related integrable system jTTj. 

The key to understand these relations should exist in a localization of the path integral 
in an effective theory on D-branes. The non-perturbative corrections to the SUSY gauge 
theory are coming from corrections of Euclidean BPS D-branes, whose world-volume does 
not contain the time direction and which are wrapping on non-trivial cycles in an internal 
space, if we embed the gauge theory into a string theoretical configuration with a suitable 
compactification. According to the effective theory on the Euclidean D-brane must 
be topologically twisted. Moreover, if we use the feature of the topological field theory, 
the partition function is reduced to a "bosonic" system after integrating out all fermionic 
fields due to the localization [T3J E] . 

In this paper, we discuss the relation between SUSY gauge theories and two-dimensional 
(bosonic) gauge theories from the first principle, that is, from the view point of the lo- 
calization. We first give a configuration of D-branes that realize an Af=2 SUSY gauge 
theory and consider the effective theory on the D-branes wrapping on the internal 2-cycle 
in section 2. We evaluate the expectation value of an operator, which corresponds to the 
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chemical potential of vortices on the cycle, and show that it reduces to the BF-theory. If 
we deform generically the D-brane effective theory by adding a "potential" term, we ob- 
tain the partition function of the so-called generalized two-dimensional Yang-Mills (gYM 2 ) 
theory [T^l EE] . This reduction is essential reason why the prepotential of the SUSY gauge 
theory and the partition function of the two-dimensional Yang-Mills theory are related 
with each other. In this section, we also discuss the relation between our analysis and 
the argument given in fTf\ IT8] where the authors derive the g-deformed two-dimensional 
Yang-Mills theory using a similar setup. 

Moreover, in section 3, we proceed the path integral of gYM 2 following the standard 
gauge theoretical method 0121]] • We will see that the system finally reduces to a "dis- 
crete matrix model" . The difference between the discrete matrix model and an ordinary 
random matrix model, which is a model for c=0 string theory and effective superpotential 
calculations by Dijkgraaf and Vafa, is that eigenvalues of the discrete matrix model are 
discretized in a unit while the eigenvalues of the ordinary matrix model are continuous 
variables. Each eigenvalue can not exist at the same position (value) due to the existence 
of the Vandermonde determinant. So a possible set of eigenvalues is given by a strongly 
decreasing (non-colliding) integer sequence. Thus the integral over possible eigenvalues 
in the partition function of the random matrix model is replaced by a summation over 
possible sets of non-colliding ordered integers. On the other hand, if we map these dis- 
crete eigenvalue distributions to a weakly decreasing sequence by a suitable shift, we can 
identify the integer sequence with numbers of rows of a Young diagram. Using this identi- 
fication, we can see the relation between the discrete matrix model and Migdal's partition 
function |21j . which is expressed as a summation over sets of the Young diagrams (rep- 
resentations of U(N)). In addition, we regard the eigenvalues as states of free fermions 
obeying the exclusion principle by using the famous Young diagram/Maya diagram (free 
fermion states) correspondence |22j . 

We find that the discrete matrix model possesses an essential structure of Nekrasov's 
instanton partition function. If we consider a potential term, the eigenvalues of the 
discrete matrix model are accumulated around the critical points of the potential, then 
each lump of the eigenvalues has two end points (fermi surfaces). We show that a suitable 
large N limit decouples these fermi surfaces with each other. Then the partition function 
is factorized into two sectors, and one of them produces the instanton partition function 
whose rank of the gauge group and moduli parameters of M=2 theory are determined by 
the critical points of the gYM 2 potential. This manner to take the large N limit is similar 
to the chiral decomposition by Gross and Taylor [213123 EE] (and see for review |26j). 

In section 4, we extend the arguments in section 3 to the M=2 theories including the 
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hypermultiplets in various representations. We first derive Nekrasov's partition function 
of theory with the hypermultiplet in the adjoint representation given in [27J E] from a 
suitable two-dimensional model. To obtain the initial two-dimensional theory, we add 
extra observables to the topological theory, and the additional terms have the same form 
as the (tree level) superpotential of the corresponding four-dimensional theory. This 
extension very looks like the Dijkgraaf-Vafa's construction to determine the matrix model 
actions, which is simply obtained from the tree level superpotential by replacing the 
superfields with the hermitian matrices. Then we get Nekrasov's partition function of this 
model from the large N limit again. Once we obtain the manner to make the observables 
for adjoint matters, it is easy to generalize it to the quiver theory 1 and the theory with 
the hypermultiplets in the fundamental representation, since these theories are essentially 
reproduced by flows from the theory with the adjoint matter. We exhibit an U(r) x 
U(s) ^-type quiver theory and an U(r) gauge theory with s flavors in the fundamental 
representation. 

In section 5, we discuss the relationship between the discrete matrix model and usual 
continuous matrix model more explicitly. There we see the discrete matrix model in the 
continuum and large N limit has two different type of third order phase transitions. One is 
called as the Douglas-Kazakov phase transition |2H| and another is the Gross- Witten phase 
transition [3U] . We give a correspondence between the limiting shape of the Young diagram 
and the eigenvalue density (free fermion states). We find that the Douglas-Kazakov 
and Gross- Witten phase transition are dual under electron/hall exchanges (exchanges of 
the eigenvalue and vacant positions). Combining these observations with the arguments 
in section 3, we also find the continuum and large N limit to derive the Nekrasov's 
instanton partition function and Dijkgraaf-Vafa's matrix model analysis should lie on 
different phases of the discrete matrix model. We expect that relations among various 
phases may connect non-perturbative dynamics of Af=2 and 7V=1 theories. 

2 Localization on D-brane and two-dimensional YM 
theory 

Four- dimensional gauge theories with eight supercharges can be realized in string theory 
in various ways like the geometric engineering [HU E2] arid Hanany- Witten type brane 
configuration j33[ EH]- Among them, we consider Type IIB superstring theory on R 1 ' 3 x 

1 The multi- instanton calculus for the quiver gauge theory is proposed by |28) from a point of view of 
the equivariant cohomology. 
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C x where M.± is a four-dimensional ALE space. If D5-branes are wrapped on 
2-cycles in M.^ a four-dimensional gauge theory will appear on the extra M 1,3 space 
on the D5-branes except for the compactified internal 2-cycles. The configuration for 
the four-dimensional gauge theory preserves 8 supercharges at least when the 2-cycle is 
CP 1 , which is T-dual to the Hanany-Witten Type IIA brane configuration of the four- 
dimensional M=2 SUSY gauge theory. So it is sufficient to consider the case of CP 1 from 
the field theoretical point of view, but we will treat genus of the 2-cycle as generic one 
throughout the paper. So we assume the structure near a single cycle is T*E G , where £ G 
is a Riemann surface with genus G, that is, the world-volume of the D5-branes is thought 
to be spanned along R 1 ' 3 x £ G . The gauge coupling of the four- dimensional system relates 
to the area A of S G by 

I A , . 

— = —. (2.1) 

9ym 4 9 s^s 

We set l s = 1 hereafter. 

In this article, we are interested in contribution of instantons to the partition func- 
tion of the four-dimensional M=2 SUSY gauge theory. To this end, we start from the 
worldvolume theory of D5-branes on R 1,3 x £ G described above. It is known that the 
theory admits noncommutative deformation on IR 1,3 and it does not affect to the final 
result of the prepotential [H3 inH] - According to the concept of the large N reduction due 
to the noncommutativity [371 EE EH] , the worldvolume theory on D5-branes is reduced 
to the two-dimensional gauge theory on a large number of Euclidean D-strings wrapping 
on E G . As discussed in [T2j, the low energy effective theory on D-branes wrapping on 
the cycle must be (partially) topologically twisted to preserve the SUSY. In this case, the 
topological theory should be localized on the Hitchin system, 

F sg = [W Z ,W Z ] , D- Z W Z = D Z W- Z = , (2.2) 

where W z and W z are 1-forms on S G associated with the normal bundle of S G . Notice 
that these fields correspond to the fluctuation of S G in AI4. However, we can simplify 
the analysis by giving a huge mass to W z [17J. Correspondingly, we freeze the degrees of 
freedom as W z = W z = and the system is localized on the flat connection F — 0. 

The twisted theory on the internal cycles enjoys the BRST symmetry and contains 
essentially a two-dimensional gauge field A = Aidx 1 (i = 1,2), its (fermionic) BRST 
partner A = \idx l and a complex scalar field $ in the adjoint representation of the gauge 
group U(N) which corresponds to the position of D-strings in the flat directions C. These 
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fields transform under the BRST symmetry as 2 



QA = A, Q\ = -d A $, Q$ = 0, (2.3) 

where d A & = c£<I> + [A 3>] • The BRST operator Q is nilpotent up to the gauge transfor- 
mation, and A, A and $ have the ghost (BRST) charges of 0,1 and 2, respectively. As 
mentioned above, we need to impose the condition F = dA + A A A = 0. To make this 
constraint and write down Lagrangian, we need additional BRST multiplets (^,rj) and 
(H, x) , which transform as 

Q$ = r), Qr) = [$,$], 
Q#=[$, X ], Qx = #- 

Using these BRST multiplets, the action of the twisted theory can be written in the BRST 
exact form as 

S top = ^Q[ TrE(A,\,<S>,$, V ,H, X ), (2.5) 



G 



where h 2 is a coupling constant which is proportional to the string coupling g s and 

5 = ^[$, $] + l X (ff - 2 * F) + d A $ A A. (2.6) 

If we write the bosonic and fermionic fields together as B = (A, H, $) and T = (A, Xi v)> 
respectively (note that $ does not have any fermionic partner), the partition function is 
given by a path integral over these fields; 

Z top = j VBVfv$e~^ Q '/eg^sC^.*) . (2.7) 

Here we notice that the action ()2.5|) can be obtained from the dimensional or large N 
reduction from the four-dimensional twisted M=2 gauge theory, where an adjoint scalar 
field and a Lagrange multiplier corresponding to the moduli and constraints of the normal 
direction to the 2-cycle are thrown away by turning on a huge mass [40J. 

From the variation of the partition function with respect to the coupling constant h 2 , 
we find the partition function ()2.7|) is independent of the coupling. This means that the 
path integral can be evaluated exactly in the WKB (weak coupling) limit. So the path 
integral dominates around the Gaussian integrals with the constraints. In particular, by 
integrating out H, this integral gives a flat connection constraint F = 0. Therefore, any 
physical observables must be evaluated around the flat connections and the path integral 
2 Q acts on X as a commutator or anti-commutator depending on whether X is bosonic or fermionic. 
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()2.7|) gives schematically the "volume" of the moduli space of the flat connections, namely 
the genus G Riemann surface with marked points. 

Next we consider BPS bound states on the D-strings. In the context of the topological 
field theory, we can introduce the chemical potentials for the BPS objects bounded on 
the D-strings by evaluating the vacuum expectation value of observables in the topolog- 
ical field theory ^7j. In general, observables in topological field theory are cohomology 
classes of the BRST operator, which are constructed from p-forms O p obeying the descent 
equations |30] . 



d A O = QO u d A O l = Q0 2) 

where, 

Oo = |Tr$ 2 , 
O x = Tr $A, 

2 = Tr (t<S>F + ±A A A) . 
Let us first consider the integral of 2-form O2] 

I 2 (£ G ) = jf Tr (i$F + ±\A\ 

and consider the vacuum expectation value, 

V&DfV$expi I Tr 



(2.1 



(2.9) 



(2.10) 



top 



AAA) 



(2.11) 



in the topological field theory discussed above. This observable corresponds to the chem- 
ical potential of vortices on the D-strings. 3 In fact, in the next section, we will show 
that 12(^0) introduces a coupling of the eigenvalues of $ with the instanton numbers 
corresponding to the maximally broken gauge symmetry U(1) N . The localization onto 
the moduli space of the flat connection due to the topological action does not affect to 
the physical part of the action j- J Sg (i&F + A A A), since the e.o.m. of the BF-theory 
also gives the flat connection constraint F = 0. So, after integrating out the bilinear of 
A, we find 1)2.11)1 exactly reduces to a partition function of the bosonic BF theory [TH] . 



top 







1 VAV^exp 


— / Tri$F 


. 9s is G 



[2.12) 



We can also consider zero-form observables in the topological field. Especially, we 
consider operators in the form of the trace of polynomials of $, which are actuarially 



3 In the context of the theory of D4-branes considered in |17|. ^(^g) relates to the chemical potential 
of DO-branes by T-duality. Hence the vortices corresponds to the intersection of D-strings. The authors 
thank to the referee of Physical Review D for pointing it out. 
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observables since $ is BRST closed itself. Recall that $ corresponds to the positions of 
D5-branes in the C-direction, which are the moduli parameters of vacua in Coulomb phase 
of the four-dimensional theory. To fix the moduli parameters in the four-dimensional 
theory, we add a generic (r + l)-the order superpotential for the adjoint scalar field 
manually, then the moduli parameters are expected to be fixed around the critical points 
of the superpotential [2|. 4 We can realize this procedure in the two-dimensional theory 
by deforming the expectation value (j2.11|) by the observable, 



r+l 

TrW(<$>) =^2^-Tr&. (2.13) 

3=1 3 

In contrast to the previous case without potential, we cannot claim that the expectation 
value of the observable deformed by the potential, -^(^g) + J E Tr W($)u> (where u is 
a volume form on Eg), is the same as the partition function of the deformed BF theory, 
since higher critical points corresponding to F ^ may contribute to the path integral. 
However, as discussed in [T^j, the contributions from the higher critical points are expo- 
nentially small. Thus, after integrating out A, we can again evaluate approximately the 
expectation value in the topological field theory as the partition function of a "physical" 
theory up to the contribution from the higher critical points; 



exp 



— [ Tr (i&F + A A A) — — / Tr W($)u 

9s JY, a 9s Js fi 



top 



— I Tr (i$F + W($)u 

s J 



;2.i4) 



VAD§ exp 

Fortunately, however, the e.o.m of the r.h.s. theory with respect to $ gives 

iF = W'($)w. (2.15) 

So we can expect that the above approximation becomes much better, if we take into 
account configurations only around the critical points W = 0, where the contribution 
from the flat connection F = is dominated. 

Hereafter we investigate the bosonic BF type theory deformed by the potential as an 
effective theory of the D-strings with vortices. This model is known as the generalized 
two-dimensional Yang-Mills theory (gYM 2 ) [T3], which is a generalization of the ordinary 
two-dimensional U(N) Yang-Mills theory. Indeed, if we choose a quadratic potential 

4 In order to recover the Af~2 supersymmetry, we reduce the superpotential to zero adiabatically after 
fixing the moduli parameters. 
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= ^$ 2 , the above partition function reduces to the usual two-dimensional Yang- 
Mills theory after integrating out $, 



ym 2 = / VAexp -tt^— [ d 2 xTiF 2 



(2.16) 



'YM2 

with an identification of <7ym 2 = 9s^2- 

Before closing this section, we need to mention on the relationship between our model 
and a system argued in [TH] - In the context of JHI, the adjoint scalar <£> comes from the 
holonomy of the gauge field at infinity of a fiber direction, which causes the periodicity of 
the field $ ~ $ + 2n. Then one needs to use a unitary measure for it and the model gets 
the g-deformation. On the other hand, in our model, we do not assume such a periodicity 
in the C-direction, which is associated with the vev for the adjoint scalar field $. Namely, 
a range of value of $ is non-compact and we do not need to use the unitary measure. In 
other words, our model is thought to be a decompacified limit of the periodic direction in 
the fiber. So we have the gYM 2 , not the g-deformed one. 



3 Instanton counting from two-dimensional gauge the- 
ory 

3.1 Migdal's partition function (Gaussian model) 

In prior to investigating the partition function of the gYM 2 , we derive Migdal's partition 
function j2I] for the ordinary two-dimensional Yang-Mills theory by using the Abelianiza- 
tion technique developed in PHI 120] • As mentioned in the previous section, the partition 
function of the two-dimensional Yang-Mills theory can be written as the gYM 2 with the 
quadratic potential, 

Z= fvAV^exp —! d 2 zTr + -$ 2 ) . (3.1) 

We now decompose the Lie algebra valued fields $ and A into 

$ = ®i T i + J2 A = J2A i T i + J2 A « E <x, (3-2) 

i a i a 

where Tj's and E a : s are the Cartan subalgebra and the root, respectively. We introduce 
the following Faddeev-Popov determinant corresponding to the gauge fixing $ a = 0; 

= Det- 1 {iJ^aM*)), (3-3) 
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where is a root vector determined by [Tj, E a ] = aiE a . Then we get 



VAfDA n V$fDc n Vc n e~~v 



pC-SGF + Sghost) 



where 



GF 



s, 



ghost 



C — (- y C/- y • 



(3.4) 

(3.5) 
(3.6) 



We first integrate out one forms A a and complex scalars (c_ a ,c a ). According to the 
Hodge decomposition theorem, any p-form £ p on a compact orientable Riemann surface 
E G can be uniquely decomposed as 



£ p = cia p _i + rf f /5 p+ i + 7, 



(3.7) 



where d) = — * d*, G Q p 1 (S G ), (3 p+ i G f2 p+1 (E G ) and 7 is a harmonic p-form. 

Applying this to A a and (c_ a , c Q ), we find that the number of modes of v4 a and (c_ a , c a ) 
is 2dimfi°(E G ) + dimi? 1 ^) and 2dimf] (E G ) + 2 dim#°(E G ), respectively. So, all 
non-harmonic modes cancel each other out and only harmonic zero modes, whose number 
is equal to 2 dimif°(E G ) — dimif 1 (E G ) = x(^ G ), contribute to the path integral in ()3.4j) . 
So we obtain the Abelian gauge theory with the partition function, 



i i^j 



Qj) 1 G exp 



(3.8) 



where we use x(£<?) = 2 — 2G. 

Recall that the BRST transformation of $ in (|2.3|) . the localization at fixed points of 
the BRST transformations, or the Gauss law constraint, tells us that the diagonal U(l) 
parts $>i(z) are holomorphic everywhere on the closed Riemann surface E G . This means 
that net effect of the path integral over is coming from an independent part in z. 
Therefore, we can put $i(z) = Aj, where Aj are ^-independent variables. So the partition 
function becomes 



Z= [HVAdXiHiK-X^exp -- f WiA^ + ^A 
J i ^ [ 9s Jt, g . I 



'CO 



(3.9) 



Thus the integration over the gauge fields becomes simply a summation over U(1) N non- 
trivial maximal torus bundles, which are classified by the first Chern class, 



dAi = 2np h pi G Z. 



(3.10) 
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Then the combination plays a role of the chemical potential for vortices on ith D- 
string wrapping on the 2-cycle as mentioned in the previous section. The integration over 
Ai reduces to the summation over a set of integers pi (i — 1, . . . , N) as 



k p k & i<j 



w - Xj) 2 ~ 2G ex.p 



1 



2%i\ k p k + 



M*2 



fe ^ i<j 



A,) 



2-2G 



- <7 s njj exp 



- n k ) exp 



n fe gZ i<j 



(3.11) 



where the factor 1/iV! corresponds to the Weyl denominator in an identification of U(N) 
U(l) x SU(N)/Z N , and we have used the Poisson resummation formula, 



^2 ex P ( -27T2— pi J = E <y(Ai/flr a - nj). 



(3.12) 



The summation over integers rij is unconstrained, but the the case of rij = rij drops 
from the above summation due to the factor ni<j( n i — n j) 2 ~ 2G if G = 0. For G > 1, 
the partition function diverges at rii = rij, but we also simply drop these singular terms 
in the summation to obtain regular results. Therefore, we can assume that n^'s are in 
a set of strongly decreasing integer sequences n\ > n 2 > ■ ■ ■ > by using the Weyl 
permutation, too. Thus we finally get a discretized version of the random matrix model 
with the quadratic potential, 5 



Z= Yl(9 s ni - g s nj) 



2-2G 



exp 



ni>H2>— >njv i<j 



E 



717 



(3.13) 



The difference from the ordinary random matrix model is that the integral over eigenvalues 
is replaced by the summation over possible sets of integer sequences. 

So far, we have ignored the normalization factor of the above path integral. In order to 
determine it, we require that a "ground state" configuration give Z = 1. The ground state 
configuration must be a densest sequence of rij and satisfy J2i n i = (setting the "origin" 
of sequence at zero) since the quadratic potential purposes to gather the eigenvalues rii 
at the origin as much as possible. The densest sequence means each difference between 



5 In [31], it is already suggested that the partition function of the ordinary two-dimensional Yang-Mills 
theory can be regarded as that of a discrete version of the random matrix model. 
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neighbors is 1, and from the condition Y2i n i = 0> we have rii = — i. Plugging back 
this configuration into the sum, the ground state contributes to the partition function as 



in* 

i<j 



si - 9s3 



,2-2G 



gstjA N(N^-l) 
e 2 12 



(3.14) 



This normalization factor itself has an important meaning since it is proportional to the 
volume of U(N) [42] . Therefore, the partition function including the normalization factor 



becomes 
Z = 



e n 

ni>n2>—>n.jv i<j 



rij — n, 



1 — 3 



2-2G 



exp 



A^(A^ 2 - 1) 
12 



(3.15) 



In order to make clear the relation to Migdal's partition function, where the partition 
function is expressed as a summation over irreducible representations of the gauge group 
U(N), we rewrite the above partition function further. The strongly decreasing sequence 
rii can be represented by a weakly decreasing integer sequence k\ > k 2 > ■ ■ ■ > through 
a relation, 

rii = ki — i + c, (3.16) 
with a constant c, since these rij's always satisfy 



rii ~ Tij = h — kj — i + j > 0, 



(3.17) 



for any i < j. Here we can regard ki as the number of boxes in i-th row of a Young 
diagram corresponding to an irreducible representation of U(N). 
Using the parametrization {ki}, the partition function reduces to 



{ki} i<j 



ki kj i + j 
i-j 



2-2G 



exp 



J2(Nki + ki(ki-2i + l)) 



(3.18) 



Notice that the dimension of the representation R associated with the Young diagram 
with {ki} is given by 

ki kj i + j 



dim /.' Yl 

l<i<j<N 



i- 3 



and the quadratic Casimir is 

C 2 (R) = J2( Nk * + Hh-2t + i)). 



(3.19) 



(3.20) 
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Figure 1: An eigenvalue distribution in the quadratic potential. Black dots stand for the 
positions of the eigenvalues and white ones are vacant positions. If we draw the right-down 
and right-up line for the black and white dots, respectively, we obtain the Young diagram 
drawing above the dots. The np and — np mean the positions of two fermi surfaces. 

Thus we finally get Migdal's partition function, 



where A = g s fiN is the 't Hooft coupling constant. Thus we now understand that the 
discrete matrix model with the quadratic potential and Migdal's partition function are re- 
lated with each other by a simple variable change, and both describe the non-perturbative 
dynamics of the two-dimensional Yang-Mills theory. 

Let us finally discuss the behavior of the eigenvalues {rii} for the gaussian model. 
We first identify the positions of the eigenvalues rij with the so-called Maya diagram. 
The correspondence between the eigenvalues of the discrete matrix model and a Young 
diagram is depicted in Fig. [T] We also regard the positions of the eigenvalues as the 
Fock space of free fermions. The states of the free fermions have two fermi surface. The 
fermions are exciting from the fermi surfaces and the corresponding Young diagram is 
scraped away (melting) from the corner. If we determine the constant c by a condition 
Ylii n i = we have c = and the fermi surfaces exist at and — The 

ground state distribution is symmetric with respect to the critical point of the quadratic 
potential, which is now at the origin. 

This choice or assumption for the ground state is also supported by the numerical 
analysis of the discrete matrix model (jH.lHj) by using Monte-Carlo simulation. For the 




(3.21) 



R 
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Figure 2: These are results from Monte-Carlo simulations for N = 100 and G — 0. 
When (3 > /3 C , the eigenvalues are accumulated around the origin exceeding the maximal 
value of the density. The corresponding Young diagram has a planar edge. If the inverse 
temperature is below the critical value (3 C , the system enters a different phase. The 
eigenvalue density behaves as the semi-circle. 

symmetric quadratic potential, the ground state is realized in the zero temperature limit 
where the potential term is dominated and the eigenvalues are symmetrically distributed. 
We can also see the fermionic eigenvalues are excited from two fermi surface if the repul- 
sive force originated in the Vandermonde determinant is taken into account. We draw 
some results of the Monte-Carlo simulation for G = in Fig|21 We vary a combination 
of the parameters (3 = g s ^A, which can be regarded as the inverse temperature of the 
system. In the low temperature (large f3) limit, we see that the eigenvalues gather at the 
origin, namely it approaches to the ground state configuration. The corresponding Young 
diagram is rigid and has a long planar edge of a rectangle shape. As the temperature 
becomes higher, the eigenvalues leave from the two fermi surfaces because of the repul- 
sive effect from the Vandermonde determinant. Correspondingly, the Young diagram is 
crumbling (melting) from the corner. On the other hand, if we define a density of the 
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eigenvalue distribution, 

1 N 

P{x) = -Y, 5 ^-^), (3-22) 
i=i 

we find that there exists a maximal value of the density reflecting the discreteness of the 
eigenvalues. A region where the eigenvalue density meets the maximal value corresponds 
to a planar (right-down) edge of the Young diagram. Looking at the behavior of the 
eigenvalue distribution (density) or the shape of the Young diagram, we notice a phase 
transition at a critical temperature. In fact, if the temperature becomes higher, the planar 
region of the Young diagram finally disappears. From the eigenvalue density point of 
view, the flat head which hits the maximal value of the density is pinched at the critical 
temperature and the so-called Wigner's semi-circle is realized in the high temperature 
phase. This type of the phase transition is known as the Douglas-Kazakov phase transition 
of third order by considering a continuum limit of Migdal's partition function or 
discrete matrix model. We will discuss it in detail in section 5. 



3.2 Nekrasov's partition function 

Here we derive the Nekrasov's instanton partition function from gYM 2 using the similar 
technique in the previous subsection. 
The partition function of the gYM 2 is 



VAV<5> exp 



-— / Tr (z$F + W($)u) 

9s JT, n 



(3.23) 



To make the eigenvalues of $ localize at suitable positions, we choose that critical points 
of the potential VT($) exist at x = qi, namely, 

r 

W'(x)=n r+1 Y[(x-q l ). (3.24) 
l=i 

The parameters q\ will be related with the moduli parameters in four-dimensional SU(r) 
Yang-Mills theory by a rescaling and we finally allow each qi to have a complex value by 
an analytic continuation. 

The gYM 2 ()3.23j) can be solved exactly by using the previous procedure for Migdal's 
partition function. The result is 

Z = J2 II te-^) 2 " 2G e-^^^-\ (3.25) 
{m} l<i<j<N 
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up to a normalization factor. So the partition function of gYM 2 can be regarded as a 
discretized version of the random matrix model with a generic potential (action). In the 
case of /i r+ i 1, the above summation over rii dominates near the critical points qi, so 
we can assume that iVj eigenvalues are around the critical point qi as 

g s ni = (<?i + g s hP, + 9sh^ ) 1 ,q 2 + g s hf ] , ...,q 2 + 9sh%], 

S v ' V v ' 

JVi N 2 

...,q r + g s h { {\ ...,q r + g s h% } r ), (3.26) 

N v ' 

N r 

where hf* (i — 1, . . . , N t ) represent fluctuations around the critical points, which satisfy 

/4°>/4°>--->4!, (3.27) 

for each /. The U(N) gauge symmetry is broken to nj=i U(Ni) around this configuration. 
In contrast to the previous Gaussian case (quadratic potential), we can not say the eigen- 
values are distributed symmetrically around each critical point in general since we can 
not ignore effects from neighbor critical points. However, in the large iV limit which we 
will discuss later, the generic potential decouples into a union of the quadratic potentials 
in an approximation. So we hereafter proceed our analysis by assuming the eigenvalues 
are symmetrical at each critical point, although this is valid only for a suitable scaling 
large N limit. 

Under this assumption, we introduce a parametrization {K^} as 

tf-Kn-i+Olll, (3.28) 

and we only consider the case 9sK® <^ \Qi — Qn\ for v n 7^ / so that the lumps of the 
eigenvalues separate enough each other. Here we have assumed iV;'s are odd numbers for 
simplicity. As well as the argument in the previous subsection, {K®} is weakly decreasing 
integer sequence; 

K[ l) > > ••• > K®, (3.29) 

and we call the densest configuration {Kf* = 0} as the "ground state" of this system. 

Let us next consider excitations from the ground state, that is, non-trivial configu- 
rations of {K^}. Since we assume symmetrical distribution of eigenvalues around each 
critical point, {K®} is assumed to be separated into two parts (positive and negative 
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Figure 3: We decompose the excitations from two fermi surfaces. We denote these exci- 
tations as kf and kf and positions of the surfaces as a\ and a\. In the large iV limit, the 
two surfaces are decoupled with each other. And one surface contributes to the instanton 
partition function. 



Kf) as 6 



K f ) = { Vm J ' , ( 3 - 3 °) 



where both of {A;^ } and } are weakly decreasing sequences of non-negative integer 

'(iV,+l)/2 



(We assume .iy 2 = ^0 ^ we cons ider the case that the numbers of the excitations 



are much smaller than Ni, most of {kf } and {A;jp} are zero and the non-zero elements 
express excitations from "fermi surfaces" at ai = qi + g s ^Y^- and on = qi — g s ^-^-. (See 
FigEl) Note that this decomposition corresponds to the chiral decomposition discussed 

in [23 121 123- 

6 As seen in (^j, this type of the limiting shape of the Young diagram is dominated in the Douglas- 
Kazakov phase (/3 > j3 c = n 2 ). 



16 



Using these notations, the potential part of the partition function becomes 



A N A r 

- W(g sni ) = - E \ W i ai + 9.(*P -i + l))+W U - g s (k? -i + 1)) 
9s i=i 9s 1=1 i=i v 7 v 7 



- E E [ w ( ai + + *)) + ^ + + *)) 9sk? 

9s 1=1 i=l 

+ W(a l - g s (-i + 1)) - W ((at - g s (-i + 1)) g s kf } ] , 

(3.31) 

where we have expanded each term by "small" fluctuations g s k\ and g s k\ . Note that 
fc^ and fc' have non-zero values only for i <C iVj, where the arguments of I4 7 ' in (J3.31|) can 
be evaluated as 

a l + « 7 ,(-i + l)~^i + 0(l). (3.32) 

Therefore, from the definition of the potential ([3.24)1 . we see that the potential part of 
the partition function is evaluated as 

r Ni 



1=1 1=1 

JVj-l 

r 2 



exp 



A 

9s 



n- 



1=1 i=\ 



C exp 



1=1 ^ ' 



(3.33) 



where Co is the contribution from the ground state; Cq = e 
and we have used a suitable analytic continuation in the anti-chiral sector. 

On the other hand, the discrete Vandermonde determinant (measure) part of the 
partition function is decomposed as 

n^.rii - g s nj) = //^V&^+V— 



(3.34) 



where 



(TV) 



MS? 



(JV) 



(AT) 



{Nl,n/i} 

IJ ( ai ~ a n + 9s{k { ^ - kfj -ii + j n )j , 
1 1 (ai ~a n + g s {k^ + kf^ - i, - j n + 1) 

(Mi)^("Jn) 

{Ni,„/2} 

IT ~ ~ 9s(Hf + - H - jn + 1) 

W,„/2} 

I I (ai - On - - ^ - U + jn) 



(3.35) 
(3.36) 
(3.37) 
(3.38) 



(J,tj)^(nj'n) 
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where the symbol {N^ n /2} over the product represents that each i\ (j n ) runs from 1 to 
( JV " 2 ~ 1 ), and we have used the suffix (N) to stress that they are defined for finite N. 
Here we find that the cross terms [A - of are decoupled from the partition function in the 
large N limit since a\ — a n ~ 0(1), ai — a n ~ 0(1) and a; — a n ~ 0(g s N). 

Now we would like to take a large iV limit. In our case, we must also take the limit of 
/i r+ i — » to recover A/"=2 SUSY in the four- dimensional space-time, thus it is proper to 
take the large N limit with iVj = iV n for all / = 1, • • • , r since the positions of the critical 
points are in the flat direction of the M=2 SUSY theory. Then, we set Ni as 

N l = 2M+l, (Z = l,---,r) (3.39) 

and take the limit M — > oo. In this limit, can be evaluated as 

(N) _ ~ / x tt / o« - «n + ^ s (4 - fc£ } - ii + jr. 



//++ = ^o(« 




0* - On + g S (~il + jn) 



with 



^o(a)= JJ (ai-a n + ^ s (-i+j)), (3.41) 

(M)^(n,j) 

where A; = Yli=i Si aric ^ * an d J now run from 1 to oo. For the proof from the first 
line to the second line of 1)3.41)1 . see the appendix A. This normalization factor is coming 
from a contribution of the ground state to the Vandermonde measure, and will give the 
perturbative part of the prepotential of the four-dimensional Af=2 SUSY gauge theory 
[HUE Elj. In addition, we take a scaling limit /x r+i — > with fixing, 

{g s M) 2r exp (-Afir+i(g.M) r ) = A 2r . (3.42) 

On the other hand, we can also show that /x gives the same contribution with 

under the setup ([3.39)1 since Hi — a n and ai — a n coincide in the large N limit ()3.39)1 . 
Therefore, in this double scaling limit, the partition function of the gYM 2 becomes 

ZgYM 2 = Z inst (a; A) 2 (3.43) 

with 

;«iA) = *(«)»» || + '' ) ia , (3.44) 
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In particular, if we set G = 0, Z- mst coincides with the Nekrasov's partition function 
exactly. 

Finally we note that, if we introduce the #-term in the action of the gYM2, we can easily 

show that the "QCD scale" A in ()3.42j) becomes complex and A appears in \i Then, 

we expect that and fi describe the instanton and the anti-instanton contributions 

to the prepotential, respectively, and we could verify 

Z 9 ym 2 = \Z imt (a; A) | 2 , (3.45) 

by taking a suitable analytic continuation of a. 



4 Adding matters 

In this section, we derive the instanton partition function for other J\f—2 theories from 
suitable two-dimensional theories. We first consider the M=2 theory with massive hyper- 
multiplets in the adjoint representation, and next we present a formula for quiver gauge 
theory [2E] as a conjecture. We also show that a formula for theory with fundamental 
hypermultiplets given in pQ is obtained by a flow from the quiver theory. 



4.1 Addition of adjoint matter 

To construct a two-dimensional theory that provides the instanton counting of the M=2 
theory with the massive adjoint hypermultiplet, we first introduce an additional BRST 
multiplets (X,i[>x) an d {Y,ip Y ) whose (deformed) BRST transformation is given by 



QX = ^j x , Qij x = [$, X]+mX, 
QY = ipy, Q^Y = [$,Y]-mY, 



(4.1) 



where X and Y are both Hermitian in the adjoint representation of U(N). Using them, 
we add the following BRST exact expression to the action of the topological field theory 

(EH), 



SS^ = / Q (Y^j x - Xt(j Y ) 

= - [ ($[X, Y] + mXY - ^y) , 

9s JT, G 

and consider the expectation value in this topological field theory, 
( exp -— / Tr ( i$F + -A A A + W($)w 

\ L 9s h G V 2 



(4.2) 



(4.3) 



top 
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where is again defined by (|3.24jl . As discussed in the section 2, it can be evaluated 

as the partition function of a corresponding two-dimensional gauge theory, 



J adj 



J VADQVXVY exp S^-j J Tr z$F + (W(&) + Y] + mXY^uj j 



(4.4) 

after integrating out A, ipx and ipY, up to contributions from higher critical points. We 
assume that these "matter" fields are localized at several points as impurities (sources) 
on the Riemann surface Eg as well as the Hitchin system in |43] . 

In the following, we show that (J4.4j) reproduces the instanton counting of the four- 
dimensional M—2 theory with the massive hypermultiplet in the adjoint representation 
given in [27J Q]. Fixing the gauge by <E> a = 0, as discussed in the section |3~T1 we obtain 
(I3.25J1 by performing the Gaussian integrals over the gauge fields A a and the ghost fields 
(c Q , c a ). In the case of ()4.4j) . we must also perform the integration over X and Y, 

n r n ^ n i 

/ n VXiPYii exp -- (m + A, - A,) X^ = R ^ ^ . _ , (4.5) 

1-+3 K &3 > l +3 J 

where we use the assumption that X and Y are delta-functionally localized at a point on 
Eg- Then, except for an overall factor, ()4.4j) can be evaluated as 

= yjr MlMg e-^^, (4 . 6) 

Now we take the double scaling limit N —>■ oo and fi r +i — y with the condition (J3.42j) . 
Repeating the discussion in the section 13 .2\ we can show that the partition function ()4.4j1 
becomes 

Z adj = Z i ti(a,m;A) 2 (4.7) 

with 



za=zs,^m)x;A" n ' " a " +9 ' (k '' i] 



1-G 



f, (l,i)+(n,j) 



ai-a n + g s (j - i) 



m + ai-a n + g s (j - i) 
m + a, - a n + g s {kl - k) ' + j - %) 



where 

a-G 



Z->(a,m)= n (a '; a " +g ^7 )r v (4.9) 
r, i! ..m + ai- a n + g s (j - i) 
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m — ► 



(a) (b) 

Figure 4: A brane configuration of J\f—2 theory with adjoint hypermultiplets. D4-branes 
are wrapping around a circle with a shift of the end points on a NS5-brane. The shift 
disappears in the massless limit and the D4-branes can be removed from the NS5-brane. 
T-duals of these brane configurations explain why the effective genus in the instanton 
partition function is increased by one. 

If we set G — 0, the chiral part Z^ t reproduces the instanton counting for the Af=2 
theory with the massive adjoint hypermultiplet. Here we have again used the fact that 
ai — a n can be replaced with ai — a n under the condition (|3.39p . 

We would like to comment on the massless limit of the above model. If we see the 
expressions (J4.8j) and (J4.9j) . we notice that the contribution from X and Y cancels one of 
the (1 — G) contributions from gYM 2 in the limit of m — > 0. This means that the genus 
G of the Riemann surface effectively shifts to G + 1. This phenomenon can be easily 
understood by using a brane configuration for the G = case. The brane configuration 
for M=2* (Donagi-Witten) [in| is shown in Fig0] where iV D4-branes are ending on a 
single NS5-brane and wrapping around a compact circle. The end points of the D4- 
branes are separated by m. Open strings between two end points have a minimal length 
and get mass proportional to m. They correspond to the hypermultiplets X and Y . In 
the massless limit, the theory enhances to J\f—4 since the two adjoint hypermultiplets are 
put into the A/"=4 vector multiplet. This means the D4-branes can be removed from the 
NS5-brane in the brane configuration and they are T-dual to parallel D3-branes. T-dual 
of the brane configuration maps the D4-branes to the D5-branes wrapping on S 2 and the 
NS5-brane to a deficit in the geometry. On the other hand, the D4-branes wrapping the 
compact circle without any shift is the T-dual of D5-branes wrapping on T 2 where the 
SUSY enhances to J\f—4. Therefore, we expect that two punctures on S 2 corresponding 
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to the end points on the NS5-brane join with each other and the topology changes from 
S 2 to T 2 in the massless limit, namely, the genus of the cycle is increased by one. This 
phenomenon should be extended to any genus G. If two punctures on the Riemann surface 
with genus G join in the massless limit, the genus is increased by one. This picture agrees 
with our derived formula. 



4.2 Quiver theory 

In this succeeding subsection, we construct a two-dimensional theory which is thought to 
provide the instanton counting of the A/"=2 A 2 -type quiver gauge theory. The generaliza- 
tion to an arbitrary v4 n -type quiver theory should be straightforward. 

The A 2 quiver theory is realized on two Riemann surfaces which intersect at a point. 
Bi-fundamental matters in the quiver gauge theory come from open strings connecting 
D-branes wrapping on the Riemann surfaces independently. So let us consider the system 
with Ni D-strings on Riemann surface Hq 1 with genus G\ and N 2 D-strings on Hg 2 
with genus G 2 . As discussed in the section the low energy effective theory of this 
system is a topologically twisted theory constructed from BRST multiplets (v4i,Ai,$i) 
and (A 2 , A2, $2), (5*1, Tjx) and ($2, 772), an d (Hi, Xi) and (H 2 , X2), which are in the adjoint 
representation of U(Ni) and U(N 2 ), respectively. The BRST transformations are given 
as well as ()2.3|) and (|2.4)L and the action is given by 

Sto P = Q(jo I Tr Nl E(A h \ h 771, Hi,xi) 

+ T2 [ TT N2 E(A 2 ,\ 2 ,^ 2 ,^ 2 ,r l2 ,H 2 ,X2))- (4.10) 

where Tr jy a denotes the trace over an N a x N a matrix, and the function S is defined like 
(|2.fij) . This theory is again independent of the coupling constants hi and h 2 . 

As in section |21 we first consider following observables of the topological field theory, 

e -j-A a) = exp j-1 Tr Na (i$ a F a + ^X a A\ a + W a {<b a )uu}j J , (a = 1, 2) 

(4.11) 

where u> a is the volume form on Sc a and we assume that the potentials W\{x) and W 2 (x) 
satisfy 

W[{x) = /i r+1 (x - q^) ■ ■ ■ (x - qW), (4.12) 
W^x) = v s+1 {x - q?) ■■■{x- q® ), (4.13) 



22 



_J_T 

respectively. If we evaluate the expectation value of e as 2 in the topological field theory, 
we obtain the product of Nekrasov's partition functions corresponding to pure M=2 SYM 
theories with the gauge groups SU{r) x SU(s) after taking a suitable large N limit. 

Furthermore, we would like to construct another "potential term" that produces the 
quiver-type interaction in the four-dimensional M=2 theory. To achieve this, we again 
introduce an additional BRST multiplets (B,iPb) and (B,ip^) which are iV\ x N 2 and 
N 2 x N\ matrices, respectively. The BRST transformations are given by 



QB = ipBi Q?Pb = $iB - £?$ 2 , 
QB = -i/> s , QtP § = m x - Q 2 B, 



(4.14) 
(4.15) 



)n i 2 1 



or, in component, 

Ni N 2 

ii=i j2=i 

Ni N 2 

{QB)i 2 i 1 — ~~ C0jj)t2*iJ iQ^B^^h ~ ^ ] B iiji (^l) jiii ~~ ^ ] (^2)iaja-^jaii; 

ii=l J2=l 

where «i, ji = 1, • • • , iVi and i 2 ? J2 = 1, • • • , ^2- We assume that these "matter" fields are 
also localized at points as impurities. Now let us define the topological field theory by 
adding the following BRST exact term, 



(4.16) 
(4.17) 



5S, 



quiv 



QTr [B^B + BiJ> s 



Ti N2 (B^b) - Ti Nl (b$ 2 B) - Tr N2 (^b) 



(4.18) 



29s js Gi us G2 

—I 

9s Js Gl us G2 

to the action (|4.10j) . The instanton partition function of the quiver gauge theory can be 
obtained from the evaluation of the expectation value, 



e 9s 



(4.19) 



top 



If we ignore contributions from higher critical points again, it reduces to the partition 
function of a two-dimensional theory, 



y qmv 



VA{DQ{DA 2 VQ 2 VBVB 



oxp<; - — 

9s J-E Gl 



1 



9s Js G2 



Tr Nl (i^Fi + Wi($i)a;i) + Tr Na \ B$xB) u x 
Tr N2 (i$ 2 F 2 + W 2 {Q 2 )uo 2 ) - Tr Nl (b<S> 2 b) lu 2 



(4.20) 
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where we have used the assumption that $1 ($ 2 ) takes non-zero value only on Ei (£2)- 
Repeating the discussion in the section the partition function ()4.20|) can be expressed 
as a summation over sets of integers {n^ } and {n^}; 



iVi 



nTT / (1) (2)\ _i -ilv"! ,^(5/)) -^2- 

11 (& n ii 9s n i 2 J e 9s 1 1 e 9s 

(4.21) 



x 

ii=l «2=1 



where A\ and A 2 are the areas of T,q x and Eg 2 > respectively. In deriving ()4.21|) . we have 
used the assumption that -B's are localized at points. 

Now let us evaluate (j4.21j) in the limit of N\ — > 00, N 2 —>■ 00 and fi r +i —> and 
u s+ i — > with a suitable scaling. Since it is almost the product of the partition function of 
gYM 2 on £ Gl and that on £ Ga , we can apply the same discussion to derive (|3.43jl . Namely, 
we can assume that there are 2iVj + 1 eigenvalues around and 2A^ + 1 eigenvalues 

(2) 

around q\ ' respectively, whose distribution can be expressed as excitations from the 
ground states, i.e. the densest and symmetrical distribution of eigenvalues around the 
critical points. When N\ and N 2 are sufficiently large, the excitations are thought to 
occur at the "fermi levels" at = + g s -^f— and = q^ — g s -%— ( a = 1, 2), 
which are labeled by sets of partitions ({A^}[ 1=1 , {k l ?}f 2=1 ) and ({fc£}F 1= i, {^}f 2 =i)- In 
the case of fi r +i and u s+ i — ► 0, it is natural to assume iVj = 2M a + 1 for each l a , 
and we can take the double scaling limit M x — > 00, M 2 - ^ 00 and /i r+ i — » and i/ a+1 — > 
with fixing, 

( 5s M0 2r exp {-A 1 ^ 1 {g s M 1 ) r ) = A?, 

(g s M 2 ) 2s exp (-A 2f i s+1 (g s M 2 ) r ) = Af. (4.22) 
(See the discussion around (|3.42jl .) In this limit, the partition function (|4.21|) becomes 

^ q uiv = [Z inst (a«;A x ) Z inst (a^;A 2 ) Z 12 (a«, a^; M x - M 2 )] 2 , (4.23) 
where Z inst (a; A) is given in (|3.43|) and 
Z 12 (a (1) ,a (2) ;M! - M 2 ) 

= E n II K ) -'?+#i-^)+J.(*f ) -*f ) -<iH] 



jgW ^< 2 > (Mi) (fc.fc) 



(4.24) 
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Here we have ignored terms that contain g s {M\ + M 2 ) since it suppresses the fluctuations 
by fci and k 2 . Therefore, the partition function (|4.20J1 is factorized into two parts in 
the double scaling limit. Especially, under the additional condition Mi = M 2 , it can be 

expressed as 

Viv = (a (1) ,a (2) ; Ai, A 2 ; d, G 2 f , (4.25) 



with 



ZZ: « (2) ; Ai, A 2 ; G u G 2 ) = (a« a (2) ) £ A^A* 2 

fc (1) ,fe (2) 



* n 

(Zl,n)^(ni j'i) 



n 

(^2-l2)^( n 2 J2) 



X 



(2) (2) . / • -x 

a}' -a> n ' + g s {j 2 - i 2 ) 
TT ± * 2 , ^ 1 ; ^, 4.26 



and 

n or - «s? + - h)) L ~ ui n (H 2) - a s + - ^) 



l-Gi -j-j. / ,,,, \ 
(Ix,ii)^(ni,j'i) (/2,«2)^(n2J2) 



II II K } - ^ + 9s (-ii + i,)) _1 , (4.27) 



X 

(Zi,ii) (^2,12) 



where fc a = f kf"\ We have again used the fact that can be replaced by as 
long as Mi = M 2 . We conjecture that the expression ()4.26|) gives the instanton contri- 
bution to the Af=2 A 2 -quiver theory with the gauge symmetry SU(r) x SU (s) with the 
bi-fundamental matters. 



4.3 Fundamental matters 

In this subsection, we derive the instanton partition function for A/"=2 theories with the 
hypermultiplets in the fundamental representation from a two-dimensional gauge theory. 
Since we already have the formula for the quiver theory in hand, it can be obtained by 
considering a limited case of the previous section, namely, the limit of A 2 — > 00. 
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Let us start with the partition function (j4.21|) . Before taking the double scaling limit 
(I4.22J) . we perform the limit of A 2 — > oo. Then the eigenvalues {flj 2 ^} cannot excite from 
the ground state since the potential becomes infinitely deep effectively, that is, we must 
take 

*£ = 0, %=0, (4.28) 

for all h and Putting Ai to A, G\ to G, {aj^ }[ 1=1 to {a;}[ =1 and {af^Yi 2 =i t° { m /}/=i> 
the factorized partition function (|4.2fi|) becomes 



7 



= ^ d (a,m;A;G) 2 (4.29) 

A2— ♦OO 



(0 /» • x 



with 

« («. mi A; G) = ZSf («. «»)£a* n l *~ °" + , ' ( * ! " " kT> + ' " !> 



ai - a n + g s (j - i) 

xnn a '-^ + f, ( -' +j) v (oo) 

(M) m a i~ m f + 9s (k\ -i+j) 

where 

(a, m) = J] (a, -a n + g s (j - i)) l ~ G ]J J[ (a, - m f + g s (-i + j))" 1 . 

(lAftnj) (i,t) (/J) 

(4.31) 

If we set G = 0, (J4.3U|) is nothing but the instanton counting for M=2 theory with massive 
fundamental hypermultiplets given in [I]. This fact supports our conjecture presented in 
the previous subsection. 



5 Continuum Limit 

5.1 Eigenvalue density vs profile function 

We first would like to mention on the relationship between the eigenvalue density and 
the profile function introduced in PQ. Both quantities are essentially connected through 
difference equations. The profile function expression of the discrete matrix model is useful 
to see a series expansion with respect to g s , which gives the graviphoton corrections to 
7V=2 SUSY gauge theory. 

We now define an eigenvalue density for ()3.25|) in the standard way, 

1 N 

P{x) = jj^25(x - g s ni), (5.1) 

4=1 
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which follows J dxp(x) = 1. The Vandermonde determinant part can be written in terms 
of the eigenvalue density as 

IJ (9 s ni - gsnj) 2 - 20 = exp I (1 - G)N 2 J dxdy log(x - y)p{x)p{y) 1 , (5.2) 

l<i<j<N I J J 

where the integral stands for a principal part. Similarly, the potential part is 

TyW($) = N J dxW{x)p{x). (5.3) 
Thus we can write the partition function of the discrete matrix model ()3.25|) as 

Z = ex p{(! - g ) n2 ~[ dxdy \og(x - y)p(x)p{y) : - / dxW(x)p(x)}, (5.4) 

{p(x)} J 9s 

where {p(x)} denotes the possible sets of the eigenvalue density p(x). 

Here we introduce a function 7 9s (x) which obeys the following difference equation with 
respect to the finite parameter g s ; 

lg s { x + 9s) + 7 Ss - 9s) - 2 79s (x) = logs. (5.5) 
Then the integral in (|5.2j) can be rewritten as 

-j- dxdy log(x - y)p(x)p(y) 

-j- dxdy(i Bm (x ~y + g s ) +lgA x - V ~ 9s) - ^lg s {x - y) S jp{x)p{y) 

dxdy<y gs (x - y)Ap(x)Ap(y), (5.6) 



where 



Ap(x)=p(x+^)-p(x-^), (5.7) 

and we have used the fact that the integral measure is invariant under a constant shift 
x — > x ± g s j1. Similarly, the potential term can be expressed as 



TrVT($) = -N j dxU{x)Ap(x), (5.1 



using the function U(x) satisfying the relation, 



U(x + ^)-U(x-^) = W(x). (5.9) 
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For r-th order polynomial potential, we generally have the (r + l)-th order polynomial 
U(x). For example, for the quadratic potential W(x) = ^x 2 , U(x) is given by 

U( x ) = £ ( ^-x 3 - % x) , (5.10) 



2 \3g s 12 

up to an integral constant. Then the partition function (|5.4jl is expressed as 
Z = exp j - (1 - G)iV 2 -j- dxdyy gs (x-y)Ap(x)Ap(y) + — J dxU(x)Ap(x)}, 

(5.11 



using the expression Ap(x). 

By the way, there is another useful instrument to describe the discrete matrix model; 
the profile function of the Young diagram utilized in pQ, 

oo 

f k (x;g s ) = \x\ + y^\\x- g s {ki -i + 1)| - \x - g s (ki - i)\ + \x + g a i\ -\x + g s (i - 1)| 
i=i 

oo 

= ^2[\x - g s (ki-i + l)\-\x - g s (ki-i)\ , (5.12) 
»=i 

where {fcj} is a sequence of non-increasing non- negative integers. We assume that hi = 
for i > N. This function is closely related to Ap(x) introduced above. To see it, let us 
consider the case where the discrete eigenvalues {n,;} are decomposed into r decoupled 
lumps starting from ai (/ = 1, • • • , r). The corresponding eigenvalue density is 

= ^ E E * (* " fl < ~ & ( k n - 1 + I) ) > ( 5 - 13 ) 

1=1 n=i ^ ^ ' ' 

where N = Ni + • — h N r , {kf^}^ l =1 is a sequence of non-increasing non-negative integer, 
and 1/2 in the delta function is a convention. Particularly, if we define the eigenvalue 
density for the "ground state" by setting all to be zero, 



i=l ii=l ^ V - / / 



we find 



Ap Q (x) = p {x + y) - p (> - y) 



1 r 

-^^(x-a^ + ^iVO -5(x-a0), (5.15) 



N 
i=i 
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since the intermediate ^-functions are canceled with each other. Using (J5.13|) . we can 
easily show that Ap(x) and the (colored) profile function is related by 

1 r 1 r 

MrS^w^ -0 ^ = A P( X ) ~ ^^2 6 ( x - a i + 9 s Ni), (5.16) 



i=i i=i 

where the second term in r.h.s. is a counterpart of the contribution to the summation 
from Ni to infinity. 

Now we have a dictionary to translate the eigenvalue density into the profile function 
in hand. As an application, let us write down the normalized Vandermonde determinant 
(measure part), 



,2 



, TT [ + - *£'+*. -<■)')_ (B . 17) 

k n-t-T -\\ ai-a n + g s {j n - i t ) > 



n\( a i9s) = exp 



using both of the eigenvalue density and the profile function. In terms of the eigenvalue 
density, (|5.17|) can be written as 

fj^(a,g s ) = exp j-iV 2 j dxdy lga {x - y)(Ap(x)Ap(y) - Ap (x)A Po (y))\ . (5.18) 

On the other hand, it can be also written as 

\ ~ \jd xd yigs{x-y)^f^ l) {x-a l )f^ n) (x~a n ) 
+ J dx^2 lg s { x ~ a i+ 9sNi) (ffan) i x - a n) -2S(x- a n )j 

l,n 

+ ^lg.(ai-0\, ( 5 - 19 ) 

using the profile function. The second term in the exponential represents the finite iV 
effect (cut-off at g 8 Ni) and the third is the perturbative contribution to the prepotential. 
These results also agree with j3] up to irrelevant constants. 

5.2 Dijkgraaf-Vafa theory and Douglas-Kazakov phase transi- 
tion 

The discrete version of the matrix model ()3.25|) has obviously a continuous limit of g s — > 0. 
In this limit, the summation over a set of possible integers reduces to an integration over 
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continuous variables if we restrict the G = (on sphere), 

r N 

z= n^n^-A^e-^^, ( 5 - 2 °) 

J i=l i<j 



where A, = g s rii and we have rescaled the coupling constant as g s /A = g s . So we have an 
ordinary type of the Hermitian (holomorphic) matrix model with the potential (action) 
W. Following [211313; we can obtain an J\f=l effective superpotential from a free energy 
of the matrix model (J5.2U)) in a large N limit by 

W eS = £ + 2mr Si, (5.21) 

i 1 i 

where T§ is a planar contribution to the free energy of the matrix model and Si = g s Ni are 
fixed for large matrix sizes iVj and identified with glueball superfields Si = Tr W^W'' 1 )" 
for i-th U(l) sector. 

This continuum limit can be understood from a point of view of the profile function 
mentioned in the previous subsection. All difference equations with the spacing g s turn 
to differential equations. For example, eqs. (|5.5|) . (|5.7|) and (|5.9|) become 

<#>J»-logz, (5.22) 
g; 1 Ap(x)^p'(x), (5.23) 
g a U'(x)^W(x), (5.24) 

in the continuum limit g s — > 0, respectively. Thus the partition function of the discrete 
matrix model in terms of the profile function ()5.11j) reduces to 

Z -> J Vp exp | (1 - G)N 2 j dxdy log(x - y)p{x)p(y) ~jj dxW(x)p(x) | , (5.25) 

where we use the relations ()5.22|) - (j5.24|) and a partial integral from the first line to the 
second line. 

Similarly, the partition functions including various type of matter discussed in the 
previous section agree with the Dijkgraaf-Vafa theory with matters in the continuum 
limit, for the adjoint matter: 

Z adj = f dMXdYe-^ [x ' Y]+mXY+wm 



TT A Jf^', Xj \ , e-^ (5.26) 



30 



the A2 quiver theory: 

= / n^II^ n ^ (A n )l = Uj) e-^ WlM+ ^ W2( ^, (5.27) 

and the fundamental matters: 



U^ WWT 1 ^ e-*K"M (5 . 28) 

i UiU f { m f + ^i) 

Thus it seems that we can reproduce all results of the effective superpotential for J\f—1 
theory from the continuum limit of the discrete matrix model, which also produces the 
prepotential for Af=2 theory with the same matter contents. 

However these continuum limits are too naive since the discrete matrix model has an 
unavoidable constraint, 

rii - n i+1 > 1, (5.29) 

for any i. Following [2E], let us consider the continuum limit g s — > and large N limit 
simultaneously, with fixing the following quantities: 

t = g s i, n(t) = g s rii. (5.30) 

The parameter t has a range of < t < S = g s N. The constraint ()5.29|) means 

n'(t) > -1, (5.31) 

in this continuum limit and it gives the constraint for the eigenvalue density: 

p(x) < 1, (5.32) 

since the n(t) and p{x) are related by 

dn _1 (a;) 1 

p{x) = — = ~wy (5 - 33) 

Therefore we must carefully take account of the constraint for the eigenvalue density when 
we solve the discrete matrix model in the continuum and large iV limit. 

This essential difference between the discrete matrix model and ordinary continuous 
matrix model causes the third order phase transition, which is called as the Douglas- 
Kazakov phase transition ■ This phase transition occurs when maximum value of the 
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eigenvalue density reaches at the upper limit 1. The Wigner's semi-circle type solution 
does not admit in this phase. 

On the other hand, there exists another type of third order phase transition in the 
ordinary continuous matrix models, namely, the Gross- Witten phase transition. Origi- 
nally, this type of the phase transition was found in a unitary matrix model, which is the 
dimensionally reduced matrix model of two-dimensional Yang-Mills theory on the lattice. 
It occurs when the end points of the lump of the eigenvalue density meet at the opposite 
side on a circle since the eigenvalues distribute on the circle due to the unit arizat ion. If 
we consider a covering space of the unitary circle, we can regard the phase transition as 
the joining of the cuts on a line of the covering space. In this sense, we can say that the 
Gross- Witten phase transition occurs when the end points of lumps of eigenvalues joins 
with each other in the continuous matrix models. 

In FigEl we sketch the shapes of the Young diagrams (the upper figures) and the 
corresponding eigenvalue densities (the lower figures) in various phases. Looking at these 
figures, we notice that the shape of each Young diagram is constructed from three pieces: 
(1) a right-down planar edge with 45 degree, (2) a right-up planar edge with 45 degree, and 
(3) a melting corner. Correspondingly, the shape of each eigenvalue density is constructed 
from (1) a "ceiling" (a flat region with the upper limit), (2) a "floor" (a flat region 
with p = 0), and (3) a curve connecting them. In fact, a right-down planar edge of 
the Young diagram corresponds to the densest eigenvalues on every g s interval, and the 
right-up planar edge of the Young diagram corresponds to the complete vacant space (no 
eigenvalues region), as typically like as the ground state distribution (a). As the corner 
of a Young diagram is melting (eigenvalues of the free fermions are exciting from the 
fermi surface), the planar edge regions get narrower. At the end, when at least one of 
the planar edge of the Young diagram disappears, a third order phase transition occurs 
as we mentioned above. The disappearance of a right-down edge means that the ceilings 
of the eigenvalue density are pinched and the Douglas-Kazakov phase transition occurs. 
Similarly, if a right-up planar edge disappears, the Gross- Witten phase transition occurs. 

From these observations, we notice that the following important fact: If we exchange 

a role of the right-down and right-up edge in the Young diagram, or correspondingly, 

the positions of eigenvalues (electron) and vacuums (hall), the Douglas-Kazakov phase 

transition and the Gross-Witten phase transition exchange each other. 7 In this sense, 

these phase transitions are essentially the same phenomena and are dual under the above 

7 For another discussion on the duality between the Douglas-Kazakov phase transition and the Gross- 
Witten phase transition, see |47|. 
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(a) Ground state 



(b) Douglas-Kazakov phase 



(c) Gross- Witten phase 



(d) Matrix model phase 



Figure 5: Various phases in the continuum limit of the discrete matrix model. Here we 
show a two-cuts solution for instance. 

operation. 

Finally, we comment on the relation to the Nekrasov's instanton partition function 
and Dijkgraaf-Vafa theory. As we discussed in section 3, we need to take a suitable 
double scaling large N limit in order to obtain the instanton partition function from the 
discrete matrix model (gYM 2 ). In this limit, the chiral decomposition is also needed. To 
decompose the partition function into two parts, we assume the melting of the corner of 
the Young diagram is sufficiently small comparing with the planar right-down edge. We 
have scaled up the planar region and decoupled the two fermi surface with each other. 
Namely, we need to take the large N limit at least in the phase (b). The scaled- up 
melting corner of one of the fermi surface will correspond to the limiting shape of the 
instanton partition function. On the other hand, the holomorphic matrix model using in 
the Dijkgraaf-Vafa theory does not have the Douglas-Kazakov type phase transition since 
there is no upper limitation in the eigenvalue densities. So all eigenvalue densities have 
the semi-circle shape. Thus, we find the Dijkgraaf-Vafa analysis is valid only for the phase 
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(d). We need to take the large N and continuum limit in this phase to get the effective 
superpotential of Af=l theories. 

We can obtain all order graviphoton corrections for the f\f—2 prepotential from an 
asymptotic expansion of Nekrasov's formula. It however is still difficult to obtain all 
order graviphoton correction for the superpotential of J\f—1 theory since we need sub- 
leading terms in the 1/N expansion of the random matrix models, where we can not use 
the WKB (saddle point) approximation. The discretization of the random matrix model 
seems to connect the two approaches by Nekrasov and Dijkgraaf-Vafa. We expect that our 
observation sheds new lights on the exact relation between them and gives a new closed 
formula for the effective superpotential of 7V=1 theory which gives all graviphoton cor- 
rections as the asymptotic expansion like Nekrasov's formula, beyond the correspondence 
between J\f—2 and M—l theories at the algebraic geometry (curve) level [IB]. 
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A Generalized ^-function and perturbative contribu- 
tion 

The perturbative contribution to the partition function is obtained essentially from infinite 
products, 

Y[ (ai ~a n + g s (j -«)), (A.l) 

(ij)GN 2 

where N is a set of natural numbers (non-zero positive integers). This type of the infinite 
product can be understood as the determinant of an operator (the product of eigenvalue 
spectra). This infinite product can be evaluated by using the (^-function regularization 
technique. 

We define the generalized ^-function as 

(a ' 2) 

where e = (ei, e 2 , . . . , ei) is a vector of parameters and n = (n 1; n 2 , . . . , no) is a set of 
non-zero positive integers. Using this ^-function, the logarithm of infinite products like 
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(|A.1|) can be expressed as 



log Y\ ( x + e 

ngN d 



n 



^2 l°g( x + e 

neN d 



Multiplying the gamma function, 



r(s) = / u^e^du, 
Jo 



by £g(s; x), we find 



du(x + e-n)-'u'- 1 e- u . 



After changing variable by it = t(x + e- ft), we obtain 



C ? (s; x) 



r(«) i, * 4 



-(x+e-n)t 



dt ts f 



-tx 



Especially, in the case of (jA.l|) . we have 

II {a l -a n + g s U-i)) = e-^- a -\ 



where 



d / i r°° dt 
—t 



-tx 



ds \T(s) J t {east _ i)( e -*t _ i) 



It is easy to find from this expression to show that it satisfies the difference equat 

l 9s ( X + 9s) + lg s { X - 9s) ~ ^lg 8 { x ) = logX, 

and that the series expansion with respect to g s becomes 

J 9s ( x ) 2 .U-r). 
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with 



fo(x) 




h{x) 



-jo logx ' 

1 1 



h{x) 



240 x 2 ' 



(A.ll) 



2g(2g - 2)x 2 9~ 2 ' 
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